GLOBAL WEAK SOLUTIONS TO THE INCOMPRESSIBLE 
NAVIER-STOKES-VLASOV EQUATIONS 



CHENG YU 

Abstract. In this paper, the system of particles coupled with fluid is considered. The 
particles are described by a Vlasov equation, and the fluid is governed by the forced 
incompressible Navier-Stokes equations. The interaction with fluid phase governed by 
Navier-Stokes equations is taken into account through a source term. The resulting 
system, namely Navier-Stokes- Vlasov equations, is shown to have global weak solutions 
in two and three spatial dimensions, and to have a unique global solution in two spatial 
dimensions. 



1. Introduction 

The fluid-particle flow has received much attention recently. It has a wide range of 
applications [5l [12l |22l [23l [25l [27j in the modeling of reaction flows of sprays, atmospheric 
pollution modeling, chemical engineering or waste water treatment, dust collecting units, 
and so on. Consequently, the description of the motion of fluid with charged particles 
has became a crucial problem. Due to the huge number of particles, it is very hard to 
calculate each particle's trajectory especially with the addition of the dimension of physical 
domain. Actually, we are not interested in the precise locations of each particle but the 
average behavior of collection of particles. For that purpose, we use kinetic method to 
describe the effect of particles on the fluid. On the other hand, macroscopic models such as 
hydrodynamic ones are commonly used in fluid. The interaction with fluid phase governed 
by fluid equations is taken into account through a force term. In this paper, we study a 
coupled system consisting of the Navier-Stokes equations of fluid dynamics and the Vlasov 
equation of the cloud of particles, namely, incompressible Navier-Stokes- Vlasov equations. 
We refer to [H H El HI [22l [231 [25l [27] for some physical introduction to Navier-Stokes- 
Vlasov equations and the related applications. 

The objective of this paper is to establish the global weak solutions for the following 
partial differential equations 

ut (u • V)u + Vp- fiAu = - (u - v)/ dv, 

(1 1) 

divu = 0, ^^■^> 
/t + v-V,./ + diVv((u-v)/) = 0, 

in X X (0, T), where C M*^ is a bounded domain, d = 2,3, u is the velocity of 
the fluid, p is the pressure, fi is the kinematic viscosity of the fluid. Without loss of 
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generality, we take /i = 1 throughout the paper. A distribution f{t, x, v) depends on the 
time t £ [0,7^], the physical position x € Q and the velocity of particle v G R'^. The 
number of particles enclosed at t > and location x € in the volume element dv is 
given by f{t, x, v) dv. 

The system is completed by the initial data 

u(0,x) =uo(x), /(0,x,v) = /o(x,v), (1.2) 

and with the following boundary conditions: 

u = on dfl, and /(t, x, v) = /(t, x, v*) for x € v • i/(x) < (1.3) 

where v* = v — 2(v • z^(x))z^(x) is the specular velocity, i^(x) is the outward normal to 0. 

When the distribution function / is absent, (jl.ip reduces to the incompressible Navier- 
Stokes equations. For the incompressible Navier-Stokes equations, we refer to the books 
[Tj El \TE[ [26] for fundamental ideas on the global weak solutions, stability, and limitations 
of the current tools. Meanwhile, there has been extensive mathematical study on the 
Vlasov equation and the related problems, for example, see [6| [T7]. 

The mathematical analysis of the coupled systems is far from being complete but re- 
cently has received much attention. Early work [16] established the global existence results 
for Stokes- Vlasov system in a bounded domain. The existence theorem for weak solutions 
has been extended in [1] , where the author did not neglect the convection term and consid- 
ered the Navier-Stokes- Vlasov equations within a periodic domain. The global existence of 
smooth solutions with small data for Navier-Stokes- Vlasov-Fokker-Planck equations was 
obtained in [13]. More Recently, the existence of global weak solutions of Navier-Stokes- 
Vlasov-Poisson system with corresponding boundary value problem was obtained in [1|. 
Meanwhile, there are considerable works in the direction of hydrodynamic limits, we re- 
fer the reader to [S] [HI [TS] [21] . In works [S] [H] [151 IH] i the authors used some scaling 
issues and convergence methods, such as compactness and relative entropy method, to 
investigate the hydrodynamic limits. A key idea in |141 [T5] is to have the dissipation of a 
certain free energy associated to the whole space, and control its dissipation rate. For the 
compressible version, local strong solutions of Euler- Vlasov equations was established in 
[2]. Global existence of weak solutions for compressible Navier-Stokes equations coupled 
to Vlasov-Fokker-Planck equations was established in [20] . 

There is no existence theory available for the Navier-Stokes- Vlasov equations in a 
bounded domain. Two different variable sets in the coupling term present new challenges 
to the existence theory. The other major difficulty arises from the Vlasov equations, such 
as the compactness properties of the distribution function /(t, x, v). The objective of this 
work is to establish the global existence of weak solutions with large data in the three 
dimensional spaces, and to establish the uniqueness in the two dimensional space. 

In this paper we shall study the initial-boundary value problem (jl.ip - (11.3p for Navier- 
Stokes- Vlasov equations and establish the global existence and uniqueness of weak solu- 
tions with large data in certain functional spaces. Motivated by the works of [HI [T9l 120] . 
we propose a new approximation scheme. The key point of approximation is to control the 
modified force term of regularized Navier-Stokes equations. The existence and uniqueness 
of the modified Vlasov equation is classically obtained, for example, see [31 [H [TB]. The 
controls of j'^a fdv and J^a v/ dv ensure that the modified force term of the Navier-Stokes 
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equations has enough regularity. Thus, we can solve the regularized Navier-Stokes equa- 
tions. The compactness properties of the system will allow us to pass the limit to recover 
the original system. We shall also establish the uniqueness of the weak solutions in the 
two dimensional space. 

The rest of the paper is organized as follows. In Section 2, we deduce a priori estimates 
for (|l.ip - (jl.3p . give the definition of weak solutions, and state our main results. In Section 
3, we construct an approximation scheme to (ll.ip - (jl.3p . establish its global existence. In 
Section 4, we establish the uniqueness for the global weak solutions in the two dimensional 
space. 



2. A Priori Estimates and Main Results 
Here we define the energy functional of the particles density: 



Hf) ■■=11 /(l+|v|2)(iv(ix. 



If u = constant, F{f) is an energy functional to the third equation in p.ip . When 
u 7^ constant, we will have the following energy inequality, more precisely: 

Lemma 2.1. The system (11. ip has an energy functional: 

E{t) := (^i|updx + F(/)) it). 

If d = 2,3 and (u, /) is a smooth solution to system (|l.ip such that 

u G L'^{0,T;L'^{n))nL'^iO,T;H^{n))- (2.1) 

/(I + |v|2) G L°^{0,T-L\n X M'^)), (2.2) 
then, for all t < T, we have: 

^E{t) = -( [ |Vupdx+ / / /|u- vpdvdx ) < 0. (2.3) 
dt \Jq JuJRd J 

Proof. Multiplying by u the both sides of the first equation in (jl.ip , and integrating over 
ft and by parts, we have 

/ -|u|2da;+ /" iVupdx = - / / /(u- v) -udvdrc. (2.4) 

Multiplying by (1 + ^|v|^) the both sides of the third equation in (|l.ip . integrating over 
and using integration by parts, one obtains that 



d_ 

dt 



F{f){t)+ I I f\u-v\'^dvdx= [ [ f{u-v) -udvdx. (2.5) 



Using (I23D-(I23]), one obtains ([23]), and ([LlD-dlS]) are the consequences of (fOjl . 
The proof is complete. 

□ 
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In what follows, we denote 



mkf= / |v|Vdv, and M^f = / |v[Vdvdx. 

Clearly, 

Mkf = / mkfdx 
Here we state the following lemmas which are due to |16] : 



Lemma 2.2. Suppose that (u, /) is a smooth solution to (jl.ip . If ^ for any p > 1, 

v)||iP < c'^^II/oIIlp, /or any t > 0; 
and if fo ^ L^{VL x W^), then we have 

1 s d+k 



|vr/dvdx<C(d,r) / IvlVodvdx +([|/o|[^.o + l)||u||^.(oT.L^+.) 

/or all a <t <T where the constant C{d,T) > depends only on d and T. 

Lemma 2.3. Under hypotheses of Lemma \2.S\ and d = 3, the density m^f and the mean 
velocity mif have the following estimates for all < t < T, 

II"^0/||l2{O) < (^(1 + ll/o||L°°{0,T;L°°{nxM3))^^ 
II"^i/||l3/2(q) < C(l + ||/o||L°°(0,T;L°°(nxR3))A^, 

where A = (J^J^s \v\^fodxdv)l + (|1/o||l°o(o,T;L°°{QxR3)) + 1) llu||i2(o,T;L6(Q))- 

Remark 2.1. Similar estimates hold in the two dimensional space. 

Our main result reads as follows. 

Theorem 2.1. For d = 2,3 there is a weak solution (u, /) to the system (jl.ip with the 
initial data (II. 2p and boundary condition (jl.3p for any T > 0. 



The weak solution to system (jl.ip - ()1.3p is defined as follows. 

Definition 2.1. A pair (u,/) is called a weak solution to the system p.ip - (|1.3p in the 
sense of distribution if 

• u G L°°(0,r;L2) nL2(o,r;i/i); 

• f{t,x,v) > 0, for any {t,x,v) e {0,T) x n x R^; 

• / G L°°(0,r;L°°(17 X R'^)nL\n X M^)); 

• /|v[2 e L°°{0,T;L^{n X M'^)); 

• for all if G C°°([0, c«) x (7) with div(p = we have 

/ / {—uipt + u-\/uip + \/u-\/ip)dxdt 
Jo Jn 

= — / f{u — v)-ipdxdvds+ / uq ■ ip{0,x) dx; 
Jo JnxR<i Jn 
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for all (f) € C°°{[0, oo) xrixM") with compact support in v, such that (f){T, x, v) = 0, 
we have 

-/ [ f{(t)t + ^-yx(t>+{^-^)-^^(t>)dxdwdt 
Jo JflxW^ 

/o0(O, X, v) dx dv. 



the energy inequality 



[l\u\'^dx+[ [ /(I + |v|2) dxdv 
Jn 2 Jn JRd. 

Vul"^ dxdt+ I I f\u-v\'^dvdxdt 
Jo JnxRd- 



< [ huo\Ux+ [ [ /o(l + |vp)dx 
Jn ^ Jn Jr'I 

holds for t e [0, T] a.e. 



In the two dimensional space, we can obtain more regularity and the uniqueness of 
global weak solution. More precisely, we have: 

Theorem 2.2. If uq G H^{n), fo S L°°{R'^)nL'^{R'^), and f^2 Wffod^^ < oo, there exists 
a unique global solution (u, /) to the system (jl.ip with the initial data p.2p and boundary 
condition (jl.3p . such that 

neL\0,T-Hlin))nC{[0,T];H^in)), — G ^^((0, T) x fi), 
/ G C{[0,T]-L°^{n X R^)) 

for any T > 0. 



3. The Existence of Weak Solutions 

The goal of this section is to show the existence of global weak solutions to p.ip with 
initial data (|1.2p and boundary condition (|1.3p . The key idea is to construct an approxi- 
mation scheme, establish its existence for the global weak solutions, and pass to the limit 
for recovering the original system. In this section, we shall prove our main result Theorem 
I2.1l in the case d = 3. All arguments do work in the case d = 2. 

3.1. Approximation Scheme. We regularize the equations (jl.ip and construct a solu- 
tion of the regularized system of equations. We view the first two equations in (jl.ip as 
Navier-Stokes equations with a source term /ig3(u — v)f dv. The key idea is to control 
/jj3(u-v)/(iv in L^{{0,T)xn) so that we can solve the Navier-Stokes equations directly. 
For that purpose, we follow the spirit in [20] to modify the Vlasov equation by truncating 
the velocity field u: we consider 

5t/ + V • VJ + divv ((xa(u) - v)/) = (3.1) 

where 

Xa(u) = Ul|lu|<;^|. 
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To preserve the similar energy inequality, we need to modify Navier-Stokes equations 
accordingly. This can be done by substituting the right hand side of the first equation in 
(ini) by 

- / (u - v)/(ivl{|u|<A}. 

To establish the global weak solutions, we find a modified Galerkin method particularly 
convenient. We define the space H as the closure of the space C^(0, M^) n {u : divu = 0} 
in L^(ri,R^). We let {4>i}'^^ be an orthogonal basis of the functional space H and such 
that 

A0j + VPi = -ei(pi in fl, 
(pi = on dn 

for i = 1,2,3.... Here < ei < 62 < 63 < ... < e„ < ... with e„ — > 00 as n — > 00. 
Define Pm '■ H H- Hm = span{(/)i, 02, ^m} as the orthogonal projection. We propose 
the following regularized Navier-Stokes equations 

dtVLm = Pm ^AUm - • Vu^ - j ^{u - ^)fm C^vl{|a|<A}^ , 

Vim{x,t) G Hm, and divu^ = 0. 

Thus the approximation scheme for the Navier-Stokes- Vlasov equations (II. ip is given by 
the following system 

dt^m — Pm (AUm ' Cm) j 

^m{x,t)eHm, and divu^, = 0. (3.2) 
dtfm + V • Vxfm. + diVv(xA(u) - v)/^) = 0, 



with the initial data, 
where 

and 



Um|t=o — -FmUo, fm\t=Q — fo, 

Gm= (U - v)/„(ivl||u|<;,| 
JR3 



u is given in L^(0,r;L2(fi)). 

The rest of this subsection is devoted to show the global existence to (j3.2p with its initial 
data. For given u, we can get enough regularity of — /]g3(u — v)/ (ivl{|u|<;^} to solve the 
modified Navier-Stokes equations. 
Indeed, we have 

Xxiu)eL^{{0,T)xn), 

due to 

u G L2(0,r;L^(O)). 
Considering the following equation 

5t/ + V • VJ + diVv((xA(u) - v)/) = 0; 

/(x, v,0) = /o(2;,v), f{t,x,v) = f{t,x,v*) ioi X e dn, V ■ u{x) < 0, 

where v* = v— 2(v-z^(x))i^(x), the existence and uniqueness of the solution can be obtained 
as in O El [16]. 
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Applying the maximal principle to the above equation, we have 

\\f{t,x,^r)\\LP <C{T)\\fo\\LP, foranyp>l. (3.3) 

Thanks to Lemma 12.21 

XA(u)G^°°((0,r)xl7), 

and 

|v|^/o dvdx < +00, 




n 

we have 

/ / \vffdvdx < +00. 
Jn Jrs 

This, together with Lemma 1 in [1], yields 

/ /dv € L2(0,r;L2(O)), and [ v/ dv € ^^(0, T; L2(17)). (3.4) 

By dSai), we get, for ah t > 0, that Gm € L'^{0, T; L'^{n)). 

For each m we define an approximate solution of (|3.2p as follows: 

m 
i=l 

and hence (|3.2p is equivalent to 

^Ml = _e,gl{t)- I {<P,{x)-Vct^u{x))Ux))dxgl{t)gi{t)+ [ G^Ux)dx. (3.5) 
The initial data becomes 



n Jn 



i=l 

which is equivalent to saying that 

5^(0) = (uo,</'i) for z = 1,2,3, ...,m. (3.6) 

So the system (|3.5|) with its initial data (|3.6|) can be viewed as an ordinary differential 
equations in verifying the conditions of the Cauchy-Lipschitz theorem. Thus it admits 
a unique maximal solution 

UmeC\[0,Tra];L\n)). 

It is easy to find the energy inequality to regularize Navier-Stokes equations as follows 

|umpc?2; + 2 / / \\7um\'^ dxdt < / \umo\'^ dx + / / Gm'^mdxdt, 
m Jo Jn Jn Jo Jn 

which, together with Gm € L^(0, T; L^(f])), allows us to take Tm = T. 
We define an operator 

S : L2((0, T) X f)) ^ L2((0, T) X 0) 

U Urn. 

Here we need to rely on the following lemma: 
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Lemma 3.1. The operator S has a fixed point in L {{0,T) x Q), that is, there is a point 
Um G L^((0,T) X il) such that Sum = u^. 



Proof. Multiplying by the both sides of (j3.2p . and using integration by parts, one 
obtains that 

^ / ^\^rn\'^dx + j |VUmpdx< j (^j ^{u - v)fm dvl{\u\<x}^ dx. (3.7) 

Considering the force term of the modified Navier-Stokes equations, we have 

1/ (u - v)/ml||u|<A}UmC?vda;| < (II / v/dv||i2 + A||/ f dv\\L2)'^ + \\u.m\\l2. 
This, together with (j3.7p . implies that 



dt 



In ^ Jn 
By Gronwall's inequality, we have 



/ ^|umpda;+ / dx < / \u.mfdx + C{m). 

Jn ^ Jn 



.up /■|u.„p..<aM 

:e(o,T) Jn 



t€{0,T) . 

which means that 

l!5'Um||L2((0,T);/fi(f7)) ^ C{m). (3.8) 

By the first equation in p.2p . one obtains that 

||5tSu^|li2(o,T;H-'{f^)) - (^-^^ 

By p.Sp and p.9p . we conclude that the operator S is compact in -L^(0, T; L^(r2)) and the 
image of the operator S is bounded in L^(0, T; fi). So Schauder's fixed point theorem will 
give us that the operator S has a fixed point in L^(0, T; i7). □ 

Applying Lemma for any T > 0, there exists a solution (um,/m) to the following 
system 



/ {Um-v)fmdvllu^l<x] 



Umix, t) £ Hm, and divu^ = 0. (3.10) 

dtfm + V • Vx/m + diVv(XA(Um) " v)/m) = 0, 

with its initial data 

Um{0,x) = PmUo, /(O, X, v) = /o(x, v), 

and boundary conditions 

Um|9Q=0, and f{t,x,v) = f(t,x,v*) for any x G v • z^(x) < 0. 
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Concerning the system (j3.10p with the initial-boundary data, we have the following 
energy inequality 

Jn ^ JnxR^ 

+ / \Vumfdxdt+ / / /„|xA(um) - vp (ivdxdt (3.11) 
Jo Jn Jo JqxR^ 

< I J|uopdx+ / /o(l + lv|2)dvdx, 
Jq ^ JnxR3 

due to the fact 

/ {Um-v)fm'i^{\u^\<x}Umdvdx = Xxi^m){Xxi^m) - fm dx dv . 

JnxR3 JnxR^ 
Then we have the following result: 

Proposition 3.1. For any T > 0, there is a weak solution (u™,/*") to the following 
system 

dtUm = Pm ^AUm - Wm ■ Vu„ - j {'^m " v)/^ (ivl|u^ | , 

u.m{x,t) G Hm, and divu^ = 0. 

dtfm + V • Vxfm + diVv(XA(Um) - v)/m) = 0, 

with its initial data 

Um(0,x) = P^uo, /(0,3;,v) = /o(x, v), 

and boundary conditions 

Umlan = 0, and f{t,x,v) = /(t,x,v*) for any x € dil, v • u{x) < 0. 
In additional, the solution satisfies the following energy inequality: 

/ l-\-a.m\'^ dx + I /m(l + [vp) dvdx 
Jn ^ JnxR^ 

/ / \Vum\'^dxdt+ / / /mlxA(um) - vp dvdxdt 
Jo Jn Jo JqxR3 



+ 



< / l\uo\^dx+ / /o(l + Iv|2)dvdx, 
Jn ^ JnxR3 



3.2. Passing to the Limit as m ^ oo. In this section, we will pass the limit as m goes 
to infinity in the family of approximate solutions (u"*,/"*) obtained in Proposition 13. li 
The estimates in Proposition 13.11 are independent of m. A, and those estimates still hold 
for any m. By (j3.3p . we have 



ll/"*llL°°(0,T;LP(nxM3)) ^ C 

for all 1 < p < cx). Proposition 13.11 yields the following estimates 

I|u'"IIl°°(0,T;L2(q)) < C, 
l|Vu'"llL2(0,T;L2(n)) < C. 
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From the above a priori estimates, we conclude that there exists a function / such that 
/"^ ^ / weak star in L°°(0, T; LP{n x M^)) 

for all p € (1, oo). 

This weak convergence cannot provide us enough information for passing the limit. For 
our purpose, we rely on the following average compactness results for the Vlasov equation 
due to DiPerna-Lions- Meyer |10] : 



Lemma 3.2. 

+ V • V^/" = diVv(F"/") in V'{rI + x (0, oo)) 



dt 

where is bounded in L°°(0, oo; P| L;^. .^,(1 + |vp)), ^"^ bounded in L°°((0,oo) x 

R^;L2(M3)j, rphf^n f^^ Pr]{v) dv is relatively compact in L'i{0,T; LP{Br)) for all R,T < 
oo, 1 < q < oo,l < p < 2 and for rj such that (^i^^^y & + L°°,a € [0, 2). 

Remark 3.1. It is crucial to use this lemma to get the strong convergence of mo/"" and 
mi/". 

Let mof and mi/ be the density and mean velocity associate with /. Applying Lemma 
T2] to the Vlasov equation of (j3.2p , one obtains that 



mof'^it, x) ^ mo/(t, x), mi/"^(t, x) ^ mi/(t, x) (3.12) 
in L'(0, T; U'{Br)) for any positive number R and 1 < q < oo, 1 < p < 2. 



Noticing that the right side of Navier-Stokes equations (13. 2p 

f (u™-v)/™dvl{|,™|<,} 

Jn 

is bounded in L°°{0,T; L^{Q,)) when A is fixed, one obtains that 

l|5tu"||L2(0,T;H-i) <C<00- (3.13) 

By (f3TT]l - (f3J3]l . applying the Aubin-Lions Lemma, (see [26]), there exist a u G L°°{0, T; L'^)f^ 
L'^{0,T;H^), such that 

u"" ^ u weak star in L°°(0, T; L^) and u*" ^ u in L^(0, T; Hh 
u™ ^ u strongly in ^^(o, T;H^). 

The next step is to show the convergence of {J^3 /"* (iv)u™'l{|um|<;^} in the sense of distri- 
butions. 

Note that Vu™ is bounded in L^{0,T, L^{n)) and (IXT2D . we have 

/"^dv) U-l{|„^|<;,} ^ (^^^ / dv^ Ul{|„|<;,| 

in the sense of distributions. Therefore, 
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in the sense of distributions. Applying these convergence results, one concludes that (u, /) 
is a weak solution to the following system 

SfU + U • Vu + Vp - Au = - / (U - v)/(ivl{[u|<;^}, 

divu^O. ^" ' 

9J + V • V,/ + diVv((XA(u) - v)/) = 

with its initial data 

m(0,2;) = uo(x), /(0,x,v) = /(0,x, v), 

and boundary conditions 

u|an = 0, and /(t, x, v) = /(t, x, v*) for x G (9il, v • i/(x) < 0. 

Next, we show that this solution satisfies a particular energy inequality. 

Because the solution (u™, satisfies the energy inequality in Proposition 13. H we have 

[ -lu™|2dx+ / /'"(! + |vl2)dv(ix 

r [ |Vu"|2dx(it+ r [ /™|xa(u'") -Vpdvdxdt 

Jo Jn Jo JnxM.3 



< [-\uo\^dx+ I /o(l + |vp)(iv(ix. 

Jo. 2 JOxRa 



The difficulty of passing the limit for the energy inequality is the convergence of the 
term /^xR^ /"^Ixa(u™) — vp dv dxdt. Here we write the term as: 



/ / /'"!xa(u™) - vpdvffxdt 
Jo JnxR^ 



r [ (nxA(u™)p-2rxA(u™)v+nv|2) dxdvdt. 

Jo Jf7xR3 



(3.16) 



By ([3T1|) . we have 

Xa(u„)^Xa(u) mL\0,T-L'^{n)). (3.17) 

Let us look at 

T 



[ [ r\xx{^n?dvdxdt- I I f\xx{u)\^ dvdxdt 
Jo JqxR'^ Jo Jqxr'^ 



<r f(f {r-f)dv)\xxiun\'ddx+r[ ([ /dv) (|xa(u™)|2-|xa(u)P) dx. 
Jo Jn \Jr^ / Jo Jn \jr3 / 

(3.18) 

Applying (f3l^ and (131711 to (I3J8I1 . we deduce that 

r [ nxA(u™)p dvdxdt ^ r I /ixa(u)p dvdxdt 

Jo JnxR'i Jo JnxR'i 

as m ^ oo. Similarly, 

[ [ vf"'x\{u"')dvdxdt^ [ [ vfxx{u)dvdxdt 
Jo JnxR^ Jo JnxRS 
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for all t > 0. 

Finally, because 

^ / weak star in L°°(0, T; LP(f] x M^)) 

for all p E (1, oo] and m2/™ is bounded in L°°[0, T; L^[U)), then for any fixed R > 0, we 
have 

r[ r\-v\^dxdwdt= r f x{\-^\ < R)\^? r dxd^rdt + 0{^) 
Jo JnxR^ Jo JnxR^ R 

uniformly in m, where x is the characteristic function of the ball of of radius R. 
Letting m ^ oo, then i? — > oo, we find 

/ [ f"'\v\'^ dxdvdt ^ [ [ f\w\^dxdwdt. 
Jo JnxR^ Jo JnxRS 

Thus, we have proved 

r[ r\xx{^n-^?dwdxdt^ r f f\xx{u)-^\^d^^dxdt (3.19) 
Jo Jf7xIR3 Jo JOxIRS 

as m ^ oo. 

Letting m go to infinity, using the convexity of the energy, the weak convergence of 
and u*", and (j3.19p . we deduce 

f-\vL\^dx+j f{l + \^r\^)dwdx 

+ [ [ |Vup(ixdt+ [ [ f\x\iu) -vl'^dvdxdt 
Jo Jn Jo JnxR^ 

< [ l\uo\^dx+ I foil + \v\^)d^rdx. 
Jn ^ JnxR3 

Thus, we have proved the following result: 

Proposition 3.2. For any T > 0, there is a weak solution (u''*,/'^) to (j3.15p with the 
initial data 

u\0,x) = uo(x), /^(0,x,v) = /o(x,v), 

and boundary condition 

u^\q^ = 0, and f^{t,x,v) = f^{t,x,v*) for x e dn, v ■ i/(x) < 0. 
In additional, the solution satisfies the following energy inequality: 

[ -\u^\'^dx+ [ f^{l+\v\'^)dvdx 
Jr^ 2 JnxR3 

+ [ [ \Vu^\'^dxdt+ [ [ /^Ixa(u^) - vpdvdxdt 
Jo Jn Jo JnxR^ 



< I J|uopdx+ / /o(l + |vp)dvdx. 
Jn 2 JnxR^ 
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3.3. Passing the limit as A — > oo. The last step of showing the global weak solution is 
to pass the limit as A goes to infinity. First, we let (u"^,/^) be a solution constructed by 
Proposition 13.21 It is easy to find that all estimates for (u™',/'") still hold for (u'^,/'^). 
So we can treat these terms as before. 

It only remains to show that we can pass the limit in the coupling terms Xa(u^)/''* and 
/r3 /'^ <^'^1{|u^|<a} = /i83 f^d'vxx{u.'^)- Here, we treat these terms as follows 

/ /V^dvl{|,.|<;,}= / /V^dv-/ /^U^dvl||,.|>;,}, (3.20) 

and for the second term in ()3.20p . 

II / f^U^ dvln^xiyx}\\L^O,T;n) 

dv||^oo(o,T;L2(n)l|u'^llL2{0,T;L6(f7))l|l{|u-^|>A}llL2{0,T;Le(n)) 

^ II f f^^..ll ll„^ll Ju^llL2(0,r;L6(n)) , 
Jr3 ^^lli°°(0,T;L2{n)l|U IIl2(o,T;L6{Q))1 j 

C 

as A — )• oo, where we used Sobolev embedding theorem. 
On the other hand, we have 



< 



(3.21) 



dt{[ f^d^r)+dw,{[ v/^(iv) = 0. 



which implies that dt{J^3 dv) is bounded in L^{0, T; H ^). This, with the help of Vu^ 
bounded in L'^{{0,T) x Q), yields that 

( / f^dv] [ fdv] as A ^ oo (3.22) 



in the sense of distributions. 

By (I3:20l) - (l3:22]) . one deduces that 



u^(/^,/^^vl„„,,„)^u(/^^y.v 



as A ^ oo 



in the sense of distributions. Thus we can pass the limit in the weak solutions of (II. ip as 
A — )• oo. We remark that the solution (u^, f^) satisfies the energy inequality in Proposition 

f-\M^\^dx+f f^{l + \w\^)dwdx 
Jn 2 JnxRS 

[ f \VM^\^dxdt+ f [ f^lxxiu^) -vl'^dvdxdt 
Jo Jn Jo JnxRS 



< [ l\uo\^dx+ [ foil + \v\^)dvdx. 
Jn ^ JnxR3 
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Using the same approach as in last subsection, letting A go to infinity, using the convexity 
of the energy and the weak convergence of and u^, we deduce 

f-\\y\^dx+f f (I + dw dx 

Jn 2 JnxR^ 

+ I I \V\i\^dxdt+ I I f\u-v\'^dvdxdt 
Jo Jn Jo JnxM.^ 

< [ l\uofdx+ [ /o(l + |v|2)dvdx. 
Jq ^ JnxR3 

So we have proved Theorem 12.11 

4. Uniqueness in The Two Dimensional Space 

The goal of this section is to establish the uniqueness of global solutions in the two 
dimensional space. For that purpose, we shall study the regularity first. 

4.1. Regularity. The existence of global weak solution to (jl.ip was obtained by Theorem 
12. li We multiply the first equation of (jl.ip by ut and use integration by parts over $7 to 
obtain. 



(4.1) 



ut\ dx + dt / |Vu| dx 
Jn 

^ / I'TT'o/ ■ u • I + / [u • Vu • Uf [ + / \mif\\ut\dx. 
Jn Jn Jn 

By the Cauchy-Schwarz inequality, we deduce 

||'5u,,2 d 2 

\\^\\LHn) + j^\\^^\\LHn) 

< ll"io/||L4(f7)||u||i4(Q)||ut||i2(f^) + ||ut||j^2(!^)||u • Vu\\l2(q) 
+ ll"il/llL2(f1)l|ut|lL2(Q)- 

By Theorem l2.lt we have 

u G L'^{0,T-L\n)),Vu € L\0,T; L'^{n)). 
Using the Gagliardo-Nirenberg inequality 

\\v\\l^ < C\\v\\%^\\Vv\\%\ (4.2) 

one obtains that 

^ /" 4 2 2 

'^^'^^ - ^II"llL°°(0,T;L2(n))li^"llL2(0,T;L2(n)) - ^- i^'^) 

Since u € L'^{0,T; Hq{Q)), by the Sobolev imbedding inequality, we obtain u € L^(0, T; L^(il)) 
for any 1 < p < oo. 

Thanks to Lemma 12.21 with d = 2, we have 

MqJ < oo for any < t < T, 

if mg/o < oo. 
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Let us estimate mo/ in the two dimensional space: 



mof= I fd^= fd^r+ / /dv 

'|v|<r •'|v|>r 



<C||/||L-r2 + l [ |v|Vrfv 

for all k >0. Letting r = {f^2 Ivj'^/fiv)^, then 

mo/<C(!|/!|L^ + l)( / Ivl'^/dv)^ 



for all A; > 0. Taking k = 6, then mo/ < C{mQfY^^ , which means 

|l"lo/||L°°(0,T;L4(n)) < oo. (4.4) 

Similarly, we have 

\Wlf\\L'-={0,T;L^(n)) < OO- (4-5) 
By (I31I-(jl3]), we have for all e > 

W^Whm + ^l|Vu||i.(^) < ^(1 + ||Vu||i.(f,)) +e||I?2u||2,(^) (4.6) 

where C{t) > 0, and C{t)dt < C for all T > 0. 

Next, observe that for all t > in view of (jl.ip . we have 

II - Au + Vp||i2(s^) 
< C (^||mi/||i2(n) + ||^||L2(f^) + |||u||Vu|||i2(f^) + ||mo/ • wWv^i^Q)^ ■ 
Due to divu = 0, by the classical regularity on Stokes equations, we obtain 

l|u|l//2(Q) 

< C (^||mi/||i2(n) + ||^||L2(f^) + |||u||Vu|||i2(Q) + ||mo/ • u||i2(Q)^ . 

Following the same argument of (|4.6p . we have for all e' > 0, 

||u||H2(n) < ^<^i(i) + C||^llL2(n) + e'\\u\\H2^n) (4.7) 

where Ci{t) > 0, and Cf{t)dt < C. Choosing e' = ^, we obtain 

||u||^2(f,)<C2(t) + C|||^||i2(f,) (4.8) 

where C2{t) > 0, and Cl{t)dt < C. Inserting (g^]) in i^M) and choose e = 1/2C, we 
obtain for all t > 

W^Whm + < C3(t)(i + l|vu||i.(^)) (4.9) 

where Cait) > 0, and Cl{t)dt < C. 

Applying Gronwall's inequality to (j4.9p . we obtain 

is bounded in L^{n x (0,T)), 
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and 

u is bounded in L°°(0,r;i7o^(Jl)). 

This, with the help ()4.8p . imphes that 

u is bounded in L2(0,r;if2(l^)). 

Here, we need to rely on the following Lemma which is a very special case of interpolation 
theorem of Lions-Magenes. We refer the readers to [26j for the proof of this lemma. 

Lemma 4.1. Let V d H <Z V' he three Hilhert spaces, V' is a dual space of V . If a 
function u belong to L'^{0,T;V) and its derivative u' belongs to L'^{0,T;V') then u is 
almost everywhere equal to a function continuous from [0, T] into H. 

Thanks to 

G L'^{n X (0,r)), and u € L°°{0,T; H^{n) D L'^{0,T- H^{n)), 
we conclude that u G C([0, T],H^{n)) by LemmaSH consequently / e CH[0, T];L°°(R'^ x 

n)). 

4.2. Uniqueness of solutions. To show the uniqueness, we rely on the following para- 
bolic regularity due to [TT| [M[26] : 

Lemma 4.2. If u solves 

dtu - Au + Vp = F, 
u{t = 0) = uo, u\qq = 0, 
on some time interval (0,T), then we have 

ll"llL°°{[0,T);Hl)nL2(0,T;H2) ^ ^ ( II -^11 L2 ({0,T) xQ) + l|uo||/fi^ 

Now we are ready to show the uniqueness. Let (ui,/i) and (u2,/2) be two different 
solutions to (ll.ip -( fL3]l . Let u = ui — U2, and / = /i — /2, then we have the following 
equations: 

Ut + Vp - Au = - (u/i + U2/ - v/) dv - (u • Vui + U2 • Vu) 

divu = 0, (4-10) 

+ V • Vj + diVv (u/i + U2/ - v/) = 

in X X (0,T), subject to the following initial data 

u(0,x) = 0, /(0,x,v) = 0, 

and boundary condition 

u|an = 0, f{t,x,v*) = f{t,x,v) for x E d^l, v • iy{x) < 0. 

Here, we denote the space X = L°°(0, T; H^) n L^{0, T; H^). Applying Lemma|M]with 
Uo = 0, we have the following regularity: 

||u|U <C\\ (u/i + U2/ - v/)(iv + (U • Vui + U2 • Vu) ||L2{(0,T)xm 
JR2 



(4.11) 



Jl + J2- 
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For Ji: 

(u/i + U2/ - v/)(iv||i2((o,T)xO) 

2 

= W^rnofi + U2mo/ - mi/||L2((o,T)xO) 

< l|u||L°°(0,T;LPi(n))ll"io/l||L2(0,T;L3{n)) + II 1 II L6 (0,T;L^(n)) II "2 ||l°° {0,r;LP2 (Q)) ||mo/ ||l3((o,T) X 
+ ll"^l/llL2{(0,T)xn) 

< sup ||mo/l||i3(n)r2 ||u||x + (^2^6 ||u2||ioo(o,T;LP(n))ll™'0/llL°°(0,T;L3(n)) 

te[o,r] 



l"^l/llL2((o,T)xn)) 



where pi = 6, x 



(4.12) 



7, for any p2 > 6, C depends on the domain And for J2: 



P2-6 ' 

||u • Vui + U2 • Vu||i2((o_r)xn) 

< lll||L'*(0,T;Ly(Q))l|u||L°°(0,T;LP{n))l|Vui||2,4((o,T)xn) 

+ l|l||L4(o,T;L°°(n)) l|u2||L4(o^T;L°=(n)) II Vu||l°° (0,T;L2(q)) 

< ||Vui||i4((o,r)xn)l|u||x + C'^*l|u2||L4(o,T;L°°(n))l|u||x, 



(4.13) 



where y = for any p > A, C only depends on the domain Q, and we used the 



Gagliardo-Nirenberg inequality for Vui, the Gagliardo-Nirenberg inequality and embed- 
ding inequality for U2. By (j4.1ip - (j4.13p . we can choose small T such that 



then 



sup ||mo/l||L3(n)T2 + C\\u2\\L°°{0,T;LPin))T(^ 
i6[0,T] 

+ CT4 ||Vui||2,4((o,T)xn) + C'^*l|u2||L4(o,r;L°°(n)) < -7^: 



|u||x < ||"1-o/||l°°(0,T;L3(J7)) + ll"T'l/llL2({0,T)xn) • 



(4.14) 



(4.15) 



The next step is to show ||w,o/||j;,3((o^T)xr2) + II™'i/IIl2((o,t)xO) can be controlled by ||u||j!s:- 
In the two dimensional space, we have 



I"^o/li3(o) dt<C 



and 







T 

Jn 

T 





|v[^/ d^r dx dt, 



|vl fd-vdxdt. 



\mlf\\L^n)dt < C 

Pluging them into (|4.15p . and choosing T small enough again, one deduces that 



n\\x<C 



Jn 



jv| fdvdxdt. 



(4.16) 
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We multiply the second equation of ()4.10p by [vl*"' for k > 1, and use integration by 
parts over Q x M^: 



k ufi\v\''-^dvdx + k / / U2f\v\''-^dvdx. 



(4.17) 



We estimate the right hand side terms of (j4.17p : 



k / / ufi\vf''dvdx + k / / U2f\vf-'dvdx 
Jn Jn JK2 



<C [ \u\\mk-ifi\dx + C [ \u2\\mk-if\dx (^•^^) 
Jn Jn 

< C'l|u||L2(n)ll"ifc-i/i||L2(n) + C'l|u2||H2(Q)||mfc_i/||ii(s^). 



By ([CTT]) and (^JEi), we have 



T 



sup / / flvfdvdx + k / / / flvfdvdxdt 
tG(o,T) Jn Jm2 Jo Jn Jm.'^ 

< \\^\\L^n)\\'mk-ifi\\L2{n)dt + ||u2||^2(s^)|[mA:_i/||Li(n)C^i 



< l|u||L2((0,T)xn)ll"T'fc~l/l|lL2(Q)T'2 + ||u2||i2(o,T;/f2(m)r2 sup / / | v]'' V dx. 

^ ° te[o,T] Jn Jr2 

(4.19) 

for all A; > 1. Meanwhile, we integrate the third equation in (|4.1U|) over nxM? and use 
integration by parts: 

/ / fdvdx= [ [ fodvdx = 0. (4.20) 
Using (|4.19p - (|4.20p and by induction, we deduce 

/ / / /Ivj'^dvdxdt < ( sup ||m3/i||i2(mri + |Iu2||x sup \\m2fi\\L2(n)T 
Jo Jn JR2 t£[o,T] te[o,T] 

+ ||u2||x sup ||mi/i||i2(mr3 + ||u2||x sup ||mo/i||L2(mr^)||u||x. 

iG[0,T] t€[0,T] 

Thanks to (|4.2ip and (|4.16p . choosing T > small enough, we obtain: 

l„-„ 

Thus, u = 0, and hence ui = U2 on a small time interval. Thus, we have proved the 
uniqueness of u on a small time interval. 

On the other hand, we have the following equation on the same time interval, 

/t + vV/ + diVv((u2-v)/) = in X X (0,r], (4.22) 

with its initial data 

/o = 0, 
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and boundary condition 

f{t, X, v) = f{t, X, V*) for X G (90, V • u{x) < 0. 

By ([02]) . we have 

II/IIl°°((0,T]xQxR2) < C'II/o||l°o((0,T]xQxR2)) 

which yields /i = /2- So we have proved the uniqueness of solution (u, /) on a small time 
interval [0,To]. For any given T > 0, we consider the maximal interval of the uniqueness, 
Ti = sup To < T, such that the solution is unique on [0,To]. The main goal is to prove 
that Ti can be taken to be equal to oo. For any given Tq > 0, we use 

u(ro,x) = 0, /(ro,x,v)=0, 

as the new initial data to (|4.10p . Applying the same argument to equation (|4.10p with the 
new data, the uniqueness of solution can be extended to [0, Tq + T*] for a small number 
T* > 0. 

By ()4.14p and (|4.2ip . T* can be chosen only depending on the upper bounds of 

sup ||mo/i||L3(n), ||ui||x, ||u2||x, sup ||mj/i||i2(Q) for z = 0, 1, 2, 3. 

tG[0,T] te[o,T] 

By the regularity of u in Section 4.1, ui, U2 are uniformly bounded in the space X. 
Applying the same argument of ()4.4p and Lemma 12.21 we can show that the other terms 
are uniformly bounded for all time t > 0. All such terms are uniformly bounded for all 
time t > 0. Thus, T* can be chosen not depending on the initial data at time Tq. In fact, 
we can choose T* = Tq. One can then repeat the argument many times and obtain the 
uniqueness of (u, /) on the whole time. 
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